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The transition is studied by means of a disorder parameter detecting condensation of magnetic monopoles in 
the vacuum. The deconfining transition is found to coincide with the chiral transition and the susceptibility p, 
related to the disorder parameter, is consistent with a first order phase transition. 



1. Introduction 

Understanding confinement and deconfinement 
in the presence of dynamical fermions is still an 
open issue in strong interaction physics. There is 
general agreement on the order-disorder nature of 
the deconfining transition in the quenched case. 
The popular order parameter is the Polyakov line 
(L); the symmetry involved is Z^. Alternatively 
the dual ('t Hooft) line (L) P can be used as a 
disorder parameter |2l3j , (order parameter of the 
disordered phase) corresponding to the dual Zn 
symmetry. 

In full QCD, i.e. in the presence of dynamical 
quarks, the situation is less clear: Zn and Zjv 
symmetries are explicitely broken. At zero quark 
mass there is a phase transition at some Tc where 
chiral symmetry is restored, the chiral condensate 
being the order parameter: at the same Tc also 
the susceptibility of the Polyakov loop shows a 
maximum which can be related to deconfinement. 

However it is not clear theoretically what the 
chiral transition has to do with the deconfinement 
transition. Moreover, in the presence of physi- 
cal quark masses also the chiral symmetry is ex- 
plicitely broken, so that neither the Polyakov loop 
nor the chiral condensate are true order parame- 
ters. 

In a series of papers j4l5l6| we have tested 
and verified in the quenched case the mechanism 
of confinement based on dual superconductivity 
of the QCD vacuum proposed by 't Hooft [7]. 
More precisely we have constructed an operator /i 



which creates a magnetic monopole in a C/(l) sub- 
group of the color gauge group selected by abelian 
projection. If the magnetic symmetry is realized 
d la Wigner and if fj, carries non zero net magnetic 
charge, then (/i) =0. Therefore (/i) ^ implies 
Higgs breaking of the magnetic U(l) symmetry 
and thus dual superconductivity. It has indeed 
been shown that the quenched QCD vacuum is 
a dual superconductor ((^) ^ 0) in the confined 
phase and goes to normal ((^) = 0) at the decon- 
finement transition. It has also been shown that 
(/i) being different or equal to zero is independent 
of the particular abelian projection chosen |6l8j . 

Both the construction of the disorder parame- 
ter (fi) and the related magnetic symmetry stay 
unaltered if dynamical fermions are introduced: 
(fj.) is therefore a good candidate disorder param- 
eter also in full QCD, especially if one expects, in 
the spirit of the Nc —^ oo limit, that the mech- 
anism which drives confinement and the decon- 
finement phase transition be the same with and 
without dynamical quarks. 

In Ref. (2 we have indeed demonstrated that 
also in full QCD with two light flavors, {^) shows 
a transition from a low temperature phase, where 
it is different from zero thus signaling dual super- 
conductivity, to a high temperature phase where 
it is exactly zero, thus signaling the disappearence 
of dual superconductivity, or deconfinement. We 
have also demonstrated that the transition for (/i) 
coincides with the chiral phase transition. Anal- 
ogous results have been obtained in Ref. by 
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measuring a parameter related to the monopole 
free energy. 

The aim of the present work is to extend the 
analysis and to show that (/i) scales with the cor- 
rect critical indices at the phase transition, thus 
demonstrating that it is indeed a valid disorder 
parameter. We will discuss the results of a finite 
size scaling analysis of (/i) around the phase tran- 
sition, showing that it gives critical indices which 
are consistent with those found through finite size 
scaling of the specific heat (and of the chiral sus- 
ceptibility), as presented in 

2. Disorder parameter 

The operator fi is defined in full QCD exactly 
in the same way as in the quenched theory |4I5I6| 



the subscript meaning the action by which the 
average is performed. In terms of p 
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Z is obtained from Z by changing the action in 
the time slice xq, S ^ S — S + AS. In the 
Abelian projected gauge the plaquettes 

IIio(^,a:^o) = 

Ui{x,xo)Uo{x + i, xo)Uj{x, xo + 6)Uq{x,xo) (2) 

are changed by substituting 

xo) ^ U^ix, xo) = U,{x, a;o)e''^^'(^^^") (3) 

where b{x — y) is the vector potential of a 
monopole configuration centered at y in the gauge 
V& = 0, and T is the diagonal gauge group gener- 
ator corresponding to the monopole species cho- 
sen. It can be shown that, as in the quenched 
case, fi adds to any configuration the monopole 
configuration b(x — y). 

Instead of (/x) we measure the quantity 



It follows from Eq. ^ that 

P=(5)s-(5)5, 



(4) 



(5) 



(A.)=expiy^ p(/3')d/3' 



(6) 



A drop of (/i) at the phase transition corresponds 
to a strong negative peak of p. 

3. Numerical Results 

We have made simulations with two degener- 
ate flavors of Kogut-Susskind quarks, using the 
standard gauge and fermion actions. Configu- 
ration updating was performed using the stan- 
dard Hybrid R algorithm. The lattice temporal 
size was fixed at Nt — 4. Different spatial sizes 
(L = 12,16,20,24,32) and values of the quark 
mass were used. For a more detailed account on 
simulation parameters we refer to |llj . 

We can assume the following general scaling 
form for (/i) around the phase transition: 

where r is the reduced temperature and m the 
quark mass. Analyticity arguments |llll2j 



sug- 
gest that in the infinite volume limit the mass 
dependence in the scaling function factorizes, so 
that P = -j^ does not depend on the mass. 

We then obtain the following scaling law: 



(7) 



In Figure 1 we show the quality of scaling as- 
suming V = 1/3, i.e. a first order phase tran- 
sition: a good agreement is clearly visible. The 
deviation from scaling in the deconfined region is 
well understood and related to the disorder 
parameter (/i) being exactly zero on that side. 

For comparison we show in Figure 2 the quality 
of scaling assuming the 0(4) critical index v = 
0.75. 

The 0(4) universality class is clearly excluded, 
while there is good agreement with a first or- 
der phase transition, confirming results obtained 
through an analysis of the specific heat and of the 
chiral susceptibility |11| . 
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Figure 1. Finite size scaling of p according to first 
order critical indices. 



Figure 2. Finite size scaling of p according to 
0(4) critical indices. 



4. Conclusions 

We have investigated the scaling properties of 
the parameter (/i), detecting dual superconduc- 
tivity of the vacuum, around the phase transition 
of full QCD with 2 flavors of staggered fcrmions, 
in order to check if (/i) behaves as a good disorder 
parameter for the confinement - deconfinement 
phase transition, as it docs in quenched QCD. 

We have shown that a finite size scaling analy- 
sis of (fi) indicates a first order phase transition 
and excludes 0(4) critical behaviour, in agree- 
ment with previous hints |9I1U| and with a de- 
tailed analysis of the specific heat and the chiral 
susceptibility jJJ. This confirms that also in full 
QCD (p) can be considered as a valid disorder 
parameter. 
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